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Abstract. The radial velocity curves of many spectroscopic binaries (SBs) are perturbed by gas streams or proxim- 
ity effects. For SBs with circular orbits, these perturbations can give rise to spurious orbital eccentricities of high 
statistical significance. But tests to identify such anomalous orbits can be constructed since perturbed velocity 
curves are in general no longer Keplerian. The derived tests are applied both to synthetic and to observed velocity 
curves. 
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1. Introduction 

Observers of spectroscopic binaries and exoplanets com- 
monly test the statistical significance of orbital eccentric- 
ities if e 0.1. The widely adopted methodology (Lucy & 
Sweeney 1971, 1973) is to test the null hypothesis Hq that 
e = 0. If Hq is rejected at the 5 per cent level, the elliptical 
orbit is accepted. If not, the circular orbit is preferred. 

Although recommending and implementing this 
methodology, Lucy & Sweeney (1971, hereafter LS) noted 
that some e's thus accepted are spurious, being due to 
absorption line distortions by gas streams or proximity 
effects. Accordingly, they suggested that, if e > is es- 
sential for a subsequent investigation, then the 1 per cent 
level of significance could be adopted. Nevertheless, they 
noted that the well known short-period eclipsing binary U 
Cephei, for which photometric evidence indicates e = 0, is 
still assigned an elliptical orbit at the 1 and even the 0.1 
per cent level of significance. 

Given that changing the level of significance in the LS- 
test is not a secure method of identifying a sample of SBs 
free from spurious e's, a different approach is required to 
detect perturbations. The expectation that a perturbed 
velocity curve will in general not be strictly Keplerian is 
the basis for the additional tests developed in this paper. 

This investigation is prompted by the work of Skuljan 
et al. (2004). Exploiting the high-precision radial veloci- 
ties attainable with a modern spectrograph, they report 
new orbital elements for the single-lined SB £ TrA de- 
rived from measurements with standard error a — 0.014 
km s _1 , a vast improvement on the previous orbit by 
Spencer Jones (1928) for which a = 1.4 km s _1 - see 
entry no. 445 in Table I of LS. With this jump in preci- 
sion, it is no great surprise that a statistically insignificant 



e of 0.06 is now replaced by a highly significant value of 
0.0140 ± 0.0002. Moreover, the other elements are corre- 
spondingly improved in precision. In particular, cr w , the 
error on the longitude of periastron u> is reduced from 63° 
to a mere 0°.9. 

With this remarkable reduction in cr w , the time inter- 
val required for the detection of apsidal motion for this 
and similar SBs is greatly reduced, raising the prospect of 
greatly increasing the data set for this classic test of stellar 
structure . But before making second epoch observations, 
an investigator should check that e, no matter how sig- 
nificant according to the LS-test, passes the further tests 
developed in this paper. 



2. Harmonic analysis 

Nowadays, orbital elements are almost exclusively deter- 
mined by least-squares fitting; and the LS-test as well as 
the tests developed herein are specific to such solutions. 
Nevertheless, these new tests are more readily understood 
in the context of determining the elements from a har- 
monic analysis of the velocity curve (Wilsing 1893, Russell 
1902, Monet 1979). 



2.1. Fourier expansion 

The Keplerian velocity V of a component of an SB with 
orbital period P is 



V — 7 + K[cos{y + u>) + e cos u>] 



(1) 
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where 7 is the radial velocity of the barycentre, K is the 
semi-amplitude, and v(t) is the true anomaly. Since the 
motion is periodic, V(t) can be expanded as a Fourier 
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V = -C Q + ^2[C 3 cosjL + Sj sinjL] 



(2) 



where L — 2n(t — Tq)/P is the mean longitude. 

For given orbital elements P, To, 7, K, ecosco, esinu), 
the Fourier coefficients Cj and Sj are derivable from 
Eq. (1) by numerical integration of the Fourier integrals. 
Alternatively, for e ^ 0.5, accurate analytic formulae can 
be obtained by expanding in powers of e. Such formulae 
are given in the Appendix for the coefficients of the first 
four harmonics. 



2.2. Estimating orbital elements 

Given numerous measurements of a component's radial 
velocity, the orbital elements may be estimated as follows. 
First, P is determined by minimizing the scatter about 
the mean velocity curve. Then, with a trial value of Tq, 
the coefficients Cj and Sj are estimated by approximating 
their Fourier integrals by summations over the observed 
velocities. The final step is to derive the orbital elements 
algebraically by equating these estimates, Cj and Sj, to 
their Kcplcrian counterparts. 

With only the leading term retained from the expres- 
sions in Appendix A, we have Co = 27, C\ = K, Si = 0, 
C'2 = Ke cosuj, and S2 = Kesinuj. Accordingly, if To 
is adjusted so that Si = 0, then 7 = Co/2, K = C\, 
e cosuj = C2/K, and esinui = S2/K. These elements can 
be improved iteratively to allow for higher order terms in 
e in the formulae for C\ t 2 and Si ,2- 

2.3. Perturbed velocity curves 

A perturbation due to gas streams or proximity effects will 
result in estimates Cj and Sj that differ systematically 
from their exact Keplerian values. But such differences in 
C1.2 and Si ,2 do not cause the above solution procedure 
to fail, and so the investigator remains unaware both of 
the perturbation and of the resulting systematic errors in 
the elements. 

To detect a perturbation, the analysis must be ex- 
tended to the third harmonic, for which the leading 
terms are C 3 = 9/8Ke 2 cos2uj and S3 = 9/8Ke 2 sin2uj. 
Evidently, C3 and S3 are predicted by the complete set 
of elements derived above from C^2 and Si^, leaving no 
adjustable parameters to fit C3 and S3. Accordingly, if the 
velocity curve is perturbed, the vector (C3, S3) will in gen- 
eral differ from the prediction (C3,Ss). This is the basis 
for the test developed in Sect. 3.2 



3. A hierarchy of tests 

In this section, tests are developed for separating SBs with 
perturbed velocity curves from those consistent with pure 
Keplerian motion. Each test compares the goodness-of-fits 



for a pair of models, where one of the pair is derivable from 
the other by imposing constraints on its parameters. 

The goodness-of-fit criterion for model k is the sum of 
weighted squared residuals 



Rh 



i(v„ - V n f 



(3) 



where v n is the observed radial velocity with relative 
weight w n , and V n is the predicted velocity. 

This sequence of tests starts with the LS-test for sta- 
tistically significant e's. For completeness and clarity, this 
test is briefly recapitulated here. 

3.1. Detection of ellipticity 

Given the data set {v n , w n ,t n ), we first calculate the least- 
squares circular orbit (k = 1), thus determining the ele- 
ments P, To, j and K. Next, the least-squares elliptical 
orbit (k = 2) is computed, thus determining the elements 
P, To, "f,K,e cosuj, esinuj. Because of the two additional 
parameters, the elliptical orbit is always a better fit to 
the data, so that R2 < R\. The question is whether the 
reduction justifies rejecting the circular orbit. To assess 
this, LS treated the problem as one of testing multivari- 
ate linear hypotheses. In this case, the null hypothesis of a 
circular orbit is obtained by adopting the following linear 
hypotheses 

e cosuj = , e sinuj = (4) 

with respect to two of the orbital elements. 

Applying standard statistical theory, LS found that the 
probability that measurement errors of an SB with e = 
could result in a least-squares elliptical orbit with R < R 2 
is 



'R 2 



(5) 



Here /?2 = 1/2 (AT — M2), where N is the number of mea- 
sured velocities, and M2 is the number of elements for the 
elliptical orbit - i.e., for the unconstrained model. 

The probability pi determines the outcome of the first 
decision we make about this SB: the circular orbit is ac- 
cepted if pi > 0.05, and the elliptical orbit if p\ < 0.05. 

If the v n were subject only to measurement errors, 
this decision would conclude the investigation. Thus, if 
Pi < 0.05, the reality of e would be accepted. But now we 
question the reality of such statistically significant e's. 

3.2. Detection of perturbation 

In Sect. 2.3, we concluded that a perturbation will in gen- 
eral cause a displacement between the vectors (C 3 , S3) and 
(C 3 , S 3 ). In the context of a least-squares solution, this im- 
plies that v n — V n , the residuals remaining after fitting an 
elliptical orbit, will still contain a 3rd harmonic compo- 
nent. 

To look for this, consider the following two-parameter 
model 

V = V e u + ACjeosjL + ASjSinjL (6) 
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Here V e u is the previously-determined elliptical orbit, and 
the additional terms are an adjustable jth harmonic, with 
j = 3 unless otherwise specified. The least-squares so- 
lution for this model (k = 3) determines the vector 
Aj = (ACj,ASj), with the best-fit subsequently written 
as 6j — (SCj,5Sj). 

Given its adjustable parameters, model 3 improves the 
fit, so that i?3 < i?2- But is this reduction large enough 
for the null hypothesis of an unperturbed elliptical to be 
rejected? Because the null hypothesis is obtained by im- 
posing two constraints, 



AC, = , AS-, 



(7) 

the required test is identical to that of LS. Thus, the prob- 
ability that measurement errors of the unperturbed ellip- 
tical orbit could give R < i? 3 for model 3 is 

with 3 = 1/2 (AT - M 3 ), where M 3 = 2 is the number of 
parameters for the unconstrained model. 

The probability p 2 determines the outcome of the first 
decision made about the reality of e. If p 2 > 0.05, the least- 
squares displacement vector S3 is consistent with measure- 
ment errors. There is then no evidence for a perturbation, 
and so e passes this first test of reality. This success implies 
that the amplitudes and phases of the first three harmon- 
ics of the velocity curve are consistent with Keplerian mo- 
tion. But if p 2 < 0.05, a perturbation is detected, and so e 
is suspect. Whether then to accept or reject e is discussed 
further in Sect. 3.4. 



3.3. Detection of Keplerian harmonic 

Fitting Eq.(6) to the velocity curve tells us that its 3rd 
harmonic can be written as if 3 = K3 + S3, where K 3 = 
(C 3 ,S 3 ). 

Independently of the significance of £3 according to 
the p2-test, we can ask if K3 is detected. This harmonic is 
regarded as detected if the hypotheses that K 3 = and 
H 3 = are both rejected. The second hypothesis is nec- 
essary in order to reject solutions that imply a fortuitous 
near cancellation of S3 and K3. 

To test the first hypothesis, K 3 is subtracted from the 
solution of Sect. 3.2 in order to see if it is essential to that 
fit's success. This subtraction is effected by again using 
Eq. (6) but now with constraints 

AC j = SC j -C j , AS, =SS j -S j (9) 

The model thus constrained is not a physically realis- 
tic model of the velocity curve. Instead, it defines a null 
hypothesis set up as a 'straw man' purely for testing the 
significance of Kj. Accordingly, if the 'straw man' sur- 
vives the test, we do not accept the null hypothesis but 
merely conclude that Kj is too small to be detected. 

With Aj fixed by Eq. (9), this model (k = 4) has no 
adjustable parameters, so that R4 is derived arithmeti- 
cally. Moreover, since Eq. (9) represents a displacement 



from the least-squares solution, R4 > i? 3 . But is this in- 
crease large enough to reject the null hypothesis? Because 
again two constraints have been imposed, the test is for- 
mally identical to that of Sect. 3.2. Thus, on the null hy- 
pothesis that model k = 4 is the true model, 



9i = 



Rj 



03 



(10) 



is the probability that measurement errors could give R < 
R3 when model 3 is fitted to the data. 

To test the second hypothesis, H 3 is subtracted from 
the least-squares solution by again using Eq. (6) but now 
with constraints 



ACj = -Cj, AS, 



Si 



(11) 



With Aj thus fixed, this model (k = 5) also has no ad- 
justable parameters, so that R§ is derived arithmetically. 
Moreover, since Eq. (11) represents a displacement from 
the least-squares solution, R$ > R3. But is this increase 
large enough to reject the hypothesis? Because again two 
constraints have been imposed, the test is again formally 
identical to that of Sect. 3.2. Thus, on the hypothesis that 
model 5 is the true model, 



, ^ 3 



03 



(12) 



is the probability that measurement errors could give R < 
R 3 when model 3 is fitted to the data. 

Given the outcomes of these two tests, we now define 



P3 = m&x(q 1 ,q 2 ) 



(13) 



which is such that p 3 < 0.05 implies qi < 0.05 and q 2 < 
0.05. Accordingly, if p 3 < 0.05, K3 is detected. But ifp 3 > 
0.05, K3 is not regarded as detected, even if q\ < 0.05. 

3.4. Decisions 

Following the pi-test's detection of a significant e, its re- 
ality is subject to the p 2 - and p 3 -tests. There are four 
possible outcomes: 

3.4.1. p 2 > 0.05; p 3 < 0.05 

The eccentricity of an SB in this domain of (p2,P3)-space 
is powerfully supported. Not only are the amplitudes and 
phases of the first three Keplerian harmonics consistent 
with the data, with no evidence of perturbation, but K 3 
is detected. Thus the interpretation of the velocity curve's 
2nd harmonic H2 as due to orbital eccentricity is de- 
cisively confirmed. Follow-up investigations that require 
e > can be carried out with complete confidence. 



3.4.2. p 2 < 0.05; p 3 < 0.05 

The eccentricity of an SB in this domain has some support. 
Despite <5 3 being significant, the detection of K 3 indicates 
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an elliptical orbit. Nevertheless, the elements are uncertain 
because of the perturbation. Follow-up investigations that 
require e > are therefore risky because of the unknown 
systematic errors of e and uj. 

3.4.3. p 2 > 0.05; p 3 > 0.05 

The eccentricity of an SB in this domain is not strongly 
supported. Although 6 3 is not significant, the failure to 
detect K 3 implies that H 2 may still be due to a 2nd 
harmonic perturbation - see Sect. 4.1. Accordingly, follow- 
up investigations that require e > are subject to risk. 

3.4.4. p 2 < 0.05; p 3 > 0.05 

The eccentricity of an SB in this domain of (p2,P3)-space 
is in doubt. The amplitudes and phases of the first three 
Keplerian harmonics are inconsistent with the data, thus 
indicating a perturbed velocity curve. Moreover, the non- 
detection of K 3 undermines the interpretation of H 2 as 
due to an elliptical orbit. Follow-up investigations that 
require e > are inadvisable. 

3.5. Comments 

The following comments are intended to clarify the above 
tests: 

a) If a circular orbit is not rejected by the pi-test, its 
acceptance is recommended. 

Some spectroscopists assume that the motivation for 
thus setting e = is a preference for simplicity; others 
believe that this decision is based only on statistical argu- 
ments. But LS (Sect. II) also stressed the physical argu- 
ment 'that the effect of tidal friction is to diminish the ec- 
centricity of a binary system.' Accordingly, when they dis- 
covered (LS, Sect. V) a high frequency of circular orbits for 
long-period SBs with giant components, they attributed 
such orbits to 'energy dissipation in the deep convective 
envelopes of the giant component.' A further discussion of 
the a posteriori evidence for the efficacy of tidal dissipa- 
tion is given by Lucy & Sweeney (1973). 

When an observer finds that e is not significant, he 
can claim, as did LS, to have discovered that a mechanism 
causing a secular decrease of e has been operating. If so, 
it is unlikely that the secular decrease has been caught 
at the special epoch when e sa E(e), the expected value 
due to measurement errors when e = - Eq. (18) in LS. 
Accordingly, if e is not significant, we expect that typically 
e <C E(e), and so the best option is to set e = 0. 

Of course, an observer may choose not to claim this 
discovery and so publish the elliptical orbit. But the re- 
sulting contamination of catalogues with e's that can be 
attributed to measurement errors is regrettable and mis- 
leading. A theorist interprets a small but non-zero e as 
indicating a modest number of e-foldings in the secular 
decrease of e. But if e is not significant, the data is con- 
sistent with an infinite number e-foldings. 



The statistical and physical arguments for setting e = 
when p\ > 0.05 were already compelling in the early 
1970's, and they are even more so now. 

b) If a circular orbit is accepted, its reality is not 
doubted. A perturbation that converts an elliptical orbit 
into an apparently circular one is possible but improbable. 

c) Because e is non-negative, errors in e are non- 
gaussian when e ~ 0. In consequence, the estimate of e 
derived from ccoslu and esinui has a positive bias, and 
this can mislead the investigator into concluding that e is 
significant. To avoid such misjudgements, LS applied the 
statistical theory of hypothesis testing. 

In this paper, the detectability of small, non-negative 
Fourier amplitudes is an issue. Again, the hypothesis- 
testing approach is preferable to comparing (biased) am- 
plitudes to their standard errors. 

d) Although pi < 0.05 is interpreted in Sect. 3.1 as 
the detection of a significant e, it could more accurately 
be described as the detection of a significant H 2 - see Sect. 
2.2. Since perturbations by gas streams or proximity ef- 
fects can introduce a 2nd harmonic even when e = 0, the 
detection of H 2 cannot of itself be decisive evidence for 
an elliptical orbit. 

e) The statistical tests of Sects. 3.1-3.3 assume 
normally-distributed measurement errors and models that 
are linear in the adjustable parameters. The required lin- 
earity holds rigorously for the p 2 - and j?3-tests - see Eq. 
(6) - but not for the pi-test. For this latter test, Monte 
Carlo experiments (Lucy 1989) confirm its validity. 

f) The adjustable 3rd harmonic in Eq. (6) is not a 
model of the perturbation. We assume only that the per- 
turbation has period P and has a Fourier expansion with 
non-zero amplitudes for one or more of the first three har- 
monics. In particular, even though we seek evidence of a 
perturbation via the third harmonic, the efficacy of the 
P2-test does not require the perturbation to contain this 
harmonic. 

To illustrate this point, let us consider an error- free 
velocity curve subject to the 1st harmonic perturbation 
eKcos L. The resulting modified elements, denoted by 
primes, arc derived using the leading-terms formulae from 
Sect. 2.2. Since only the cosine term of the 1st harmonic 
is perturbed, we have C\ = (1 + e)C\ = (1 + e)K, with 
all other estimated coefficients having their exact values. 
Accordingly, we find 7' = 7, K' = K(l + e), uu' = oj, 
and e' = e/(l + e). On the basis of these elements, we 
predict C' 3 = 6*3/(1 + e) and similarly for S 3 . Thus, as 
anticipated in Sect. 2.3, there is a discrepancy between 
the predicted vector (C 3 ,S 3 ) = (C 3 ,S3)/(1 + e) and its 
estimate (C 3 , S 3 ) = (C 3 , S3). 

g) When the p2-test detects a perturbation, the orbital 
elements are subject to unkown systematic errors. No sta- 
tistical procedure can then extract unbiased elements from 
the velocity curve. Additional data or a physical model of 
the perturbation is required to derive more reliable ele- 
ments. 
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h) Eq. (13) in LS defines a quantity 

" = iS ,14) 

which is a dimensionless measure of a velocity curve's pre- 
cision. In terms of fi, e should be detectable at the 5 per 
cent level if 

e > 2.45^ (15) 

- see Eq. (22) in LS. 

An essentially identical analysis allows criteria to 
be derived for the detection of Keplerian harmonics in 
perturbation- free orbits with e 2 <C 1. Thus, the 5 per cent 
criteria are 

e 2 > 2.18/u (16) 

for K 3 , and 

e 3 > 1.84/x (17) 

for K 4 . Generally, Kj should be detectable at the lOOp 
per cent level of significance if 

j-^ y^Bw (18) 

where aj is given by Eq. (A. 11). 

The inequality given by Eq. (16) defines the domain 
in (e, ^i)-space within which e can be confirmed by the po- 
test. Ideally, an observing campaign should be designed to 
obtain a velocity curve that falls in this domain. 

i) The p2-test assumes that the perturbation has pe- 
riod P. But some causes of perturbation violate this as- 
sumption. One example is pulsations, another is spots, 
and yet another is a third-body perturbation. In general, 
such effects increase the residuals and hence the deduced 
measurement error a and are not detected by p 2 - 

j) The ;>2-test relies entirely on radial velocities for 
detecting a perturbation. Phase-dependent line strengths 
and line shapes may provide independent evidence of per- 
turbation. Such evidence may also resolve the ambiguity 
noted in Sect. 3.4.3. 

k) Following common practice, we have adopted the 
5 per cent level of significance in stating the decisions 
reached from the values of Pi,2,3- But the level can be 
varied. In particular, the impact of erroneous conclusions 
on subsequent investigations should be considered. 

1) Many astronomers are now familiar with defin- 
ing confidence intervals about a minimum-x 2 solution in 
terms of an increment A% 2 - e.g. Press et al. (1992). The 
relevant statistical theory is related to that from which 
the tests of Sect. 3 are derived. 

In the limit N — > oo, the tests of Sect. 3 simplify to 

Ax 2 = -2£np (19) 

so that p < 0.05 corresponds to A% 2 > 5.99. Eq. (19) 
reproduces the column for v = 2 in the table on p. 692 of 
Press et al. (1992), except that their p is the enclosed, and 
ours the excluded, probability. 

Note that, in contrast to Eq(19), the tests of Sect. 3 
do not require N — * oo. Monte Carlo simulations confirm 
that they remain valid even when N only slightly exceeds 
M. 



4. Monte Carlo experiments 

In this section, the tests of Sect. 3 are applied to synthetic 
velocity curves. With the exact elements known, these ex- 
periments investigate the reliability of the statistical tests. 

4.1. Perturbed circular orbit 

Probably the commonest error for SBs is the assignment 
of a significant e to a binary with vanishingly small e but 
subject to a 2nd harmonic perturbation. To simulate this 
case, synthetic data sets (v n ,t n ) are generated from the 
formula 

v n = V n + eKcos(2L n - zu) + uz n (20) 

Here the first term V n is the exact radial velocity at t n for 
a circular orbit of semi-amplitude K; the second term is 
a 2nd harmonic perturbation of dimensionless amplitude 
e; and the third term is a random gaussian measurement 
error with variance cr 2 . 

A least-squares elliptical orbit for such a velocity curve 
will have a spurious ewe with u> w w. By repeated sam- 
pling, we can investigate the success of the pi- and p 3 -tests 
in detecting these spurious e's despite their possible sig- 
nificance according to the pi-test. 

The first experiment is as follows. For fixed e, we gener- 
ate and analyse 10 3 synthetic velocity curves with random 
w's and uniformly-distributed L n . Each data set com- 
prises N = 50 velocities v n with a/K = 0.02, so that 
/i = 4x 10~ 3 . Finally, the least-squares orbits assume P 
fixed at its exact value. From the 10 3 analyses, the frac- 
tions /i,2,3 with pi,2,3 < 0.05 are computed, and this is 
then repeated with e increasing from to 0.15. 

The results are plotted in Fig.l. The f\ -curve shows 
that significant e's increase sharply from w 5 per cent at 
e = to 100 per cent for e k, 0.02. The 50 per cent de- 
tection rate occurs close to the value ei = 2A5/J, — 0.0098 
expected heuristically from Eq. (15) - the point Pi. The 
/2-curve shows that for e Ss 0.06 only a small fraction of 
the significant e's are detected as spurious by the p2-test. 
But with further increase of e, the absence of K a in Eq. 
(20) becomes detectable, and the / 2 -curve rises to give 
100 per cent detections of £3 when e <; 0.12. The 50 per 
cent rate occurs close to the value €2 = 0.093 similarly 
expected from Eq. (16) - the point P 2 . 

Finally, the /3-curve shows that for e *J 0.09, a small 
fraction (w 0.06) of spurious e's are adjudged real by the 
P3-test. 

This experiment demonstrates that, for perturbed cir- 
cular orbits, the tests of Sect. 3 perform as expected. But 
the results also emphasize that if a significant e is such 
that 

2.45^ £ e £ ^2.18^ (21) 

then we cannot decide if e is real or spurious. In this inter- 
val, the p2-test cannot detect a spurious e, and the p3-test 
fails to confirm a real e. 
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Fig. 1. Monte Carlo simulations. The fractions /i,2,3 of syn- 
thetic velocity curves with pi, 2, 3 < 0.05 are plotted against e, 
the fractional amplitude of the 2nd harmonic perturbation in 
Eq.(20). The starred points indicate the theoretical values of e 
at which typically the spurious e becomes significant (Pi) and 
the perturbation is detected {P2). 



4.2. Perturbed elliptical orbit 

In a second experiment, e = in Eq. (20) is replaced by 
e = 0.12 with random u. Other details are unchanged. 

From a sequence of simulations with e increasing from 
to 0.15, we find that the £»2-test detects the perturbation 
with a reliability of 50 and 90 per cent at e w 0.035 and 
0.06, respectively. The p3-test confirms the reality of e 
with a reliabilty that drops from ps 95 per cent at e = 
to w 60 per cent at e = 0.15. This drop when e w e is 
not surprising since the imposed perturbation in the 2nd 
harmonic then has an amplitude w \K?\. 

4.3. Unperturbed elliptical orbits 

In a third experiment, the perturbation amplitude e = 
and e is now the varied parameter. 

From a sequence of simulations with e increasing from 
to 0.13, we find, as expected from Sect. 3.5 i), that the 
Pi-test has a <; 50 per cent succes rate in detecting e when 
e ^ 2.45/1, and the p3-test a ^ 50 per cent succes rate in 
confirming the reality of e when e 2 <; 2.18/x. 



5. Examples 

In this section, the tests developed in Sect. 3 are illustrated 
by applying them to SBs with published orbits. 



5.1. HD 45088 

This SB is the subject of the first paper (Griffin & 
Emerson 1975) in Griffin's series on spectroscopic orbits 
from photoelectric radial velocities. The authors report a 
high quality orbit with e = 0.150 ± 0.004 derived from 
N = 54 measured velocities with standard error u = 1.0 
km s _1 . The eccentricity is clearly highly significant but, 
given the relatively short period of 6"?992, one might sus- 
pect that the orbit is circular and e spurious. 

The first step is to apply the LS-test for the significance 
of e. This and subsequent calculations are carried out with 
P fixed at the value determined by the original authors, so 
that M 2 = 5. The result of the LS-test is p x = 5.6 x 10~ 38 , 
confirming the overwhelming significance of e. 

The second step is to compare this elliptical orbit to 
the fit obtained with Eq. (6) when j = 3. The least-squares 
fit gives SC 3 = 0.28 ± 0.20 km s" 1 and 6S3 = 0.07 ± 0.19 
km s _1 , suggesting little or no evidence for 63 beyond 
what can be attributed to measurement errors. This is 
confirmed by Eq. (8), which gives p 2 = 0.35. The am- 
plitudes and phases of the first three harmonics of the 
velocity curve of HD 45088 are therefore consistent with 
unperturbed Keplerian motion. 

The elliptical orbit of HD 45088 is now subject to the 
test of Sect. 3.3. The 3rd harmonic coefficients computed 
from the least-squares elliptical orbit are C3 = —1.26 km 
s _1 and S3 = 0.58 km s _1 . Subtracting this component 
according to Eqs. (6) and (9), then computing R4, we 
find that q\ = 1.2 x 10~ 7 . Then, similarly subtracting 
H 3 according to Eqs. (6) and (11), we find that q 2 = 
3.9 x 10~ 6 . Accordingly, p 3 = q 2 = 3.9 x 10~ 6 and so there 
is a highly significant detection of K 3 . This confirmation 
of ellipticity and the earlier non-detection of perturbation 
place HD 45088 in the domain of (p2,P3)-sp&ce containing 
orbits with the most reliable e's. 

Having detected K 3 , we can extend the search to K4. 
However, the predicted coefficients C\ — —0.15 km s^ 1 
and 5*4 = —0.19 km s _1 are somewhat too small for de- 
tection. 

The calculations for the 3rd and 4th harmonics are 
illustrated in Fig. 2. For the 3rd harmonic, we see that 
the components of £3 = P3Q3 are comparable with their 
standard errors, in agreement with their non-significance 
according to pi ■ It is also evident that K 3 = OP 3 is the 
dominant contribution to OQ3, so that its subtraction 
significantly increases the residuals, leading to its highly 
significant detection (p 3 <C 0.05). 

The corresponding vectors and error bars for the 4th 
harmonic are also plotted in Fig. 2. These show that de- 
tection of this harmonic is just beyond reach. 

The outcome of this investigation of HD 45088 is that 
it passes all tests that can be usefully carried out given the 
precision of the Griffin-Emerson velocity curve. Needless 
to say, new data with fi reduced significantly below the 
Griffin-Emerson value of 3.7 x 10~ 3 might reveal a pertur- 
bation that is currently undetectable. On the other hand, 



Lucy: Elliptical orbits 



7 



HD 45088 



liar m o n i c; 
coefficients 




-1.5 -1 -0.5 

Cj (km/s) 

Fig. 2. Fourier coefficients for HD 45088. The 3rd and 4th 
harmonics of the least-squares elliptical orbit are shown as the 
vectors OP3 and OP4, respectively. The displacement vectors 
63,4 obtained by fitting Eq.(6) are also plotted together with 
their la error bars. 



Harm on i c 
c oefficienLs 




< TrA 



-4 -2 
C 3 (m/s) 



Fig. 3. Fourier coefficients for £ TrA. The 3rd harmonic of 
the least-squares elliptical orbit is shown as the vector OP 3 . 
The displacement vector £3 obtained by fitting Eq.(6) is also 
plotted as P3Q3 together with la error bars. 



if a perturbation were still not found, the confirmation of 
Keplerian motion could be extended to the 4th and higher 
harmonics. 

For the above value of fi, Eq. (15) predicts that e > 
0.009 is required for the detection of e. This is consistent 
with e = 0.15 being overwhelmingly significant according 
to the pi-test. Similarly, Eq. (16) predicts that e > 0.09 
is required for the detection of K 3 , consistent with its 
highly significant detection according the p2-test. Finally, 
Eq. (17) predicts that e > 0.19 is required for the detection 
of K4, consistent with its non-detection. 

5.2. ( Trianguli Austral is 

For this 12<?976 SB, Skuljan et al. (2004) report 225 radial 
velocities obtained with a high-resolution fibre-fed spec- 
trograph. As noted in Sect. 1, the gain in precision over 
the earlier orbit by Spencer Jones (1928) is a remarkable 
two orders of magnitude. 

Recomputing the orbit with P fixed (their Table 6), we 
find p\ = 1.5 x 10~ 148 , in agreement with their claim of 
high statistical significance despite the small eccentricity 
of 0.014. 

Let us now investigate the reality of e by applying the 
tests of Sect. 3.2 and 3.3. The first step is to fit Eq. (6) to 
their data. Switching the unit of velocity to metres per s, 
we obtain 8C3 = —3.3 ± 1.3 m s _1 and 6S3 = —3.2 ±1.5 
m s _1 , which suggest a significant perturbation. This is 
confirmed by finding P2 = 0.0063. Accordingly, the null 
hypothesis of an unperturbed elliptical orbit is rejected. 

The next step is to seek confirmation of K3, for which 
the predicted components are C3 = — 1.34 m s _1 and 



S3 = 0.96 m s^ 1 . Since these are comparable to the stan- 
dard errors on the components of 63, we immediately sus- 
pect that there is little support for the detection of this 
harmonic. This is confirmed by the p 3 -test. Subtracting 
-K3, we find that q\ = 0.47. Similarly, subtracting if 3 , we 
find that q 2 = 9.9 x 10~ 4 . Accordingly, p 3 = q 1 = 0.47, 
and so K3 is not detected. 

These results place ( TrA in the domain of (p2,P3)- 
space discussed in Sect. 3.4.4. This contains the least reli- 
able category of e's. The dramatically improved precision 
has allowed the discovery of a very small amplitude yet 
highly significant if 2 But its interpretion as due to or- 
bital eccentricity is put in doubt by the detection of a sig- 
nificant S3 and the non-detection of K3. The amplitude 
of H2 ~ 105 m s _1 is so small that there is no difficulty 
in attributing it entirely to the detected perturbation. 

The above calculations are illustrated in Fig. 3. 
Comparison with Fig. 2 for HD 45088 is illuminating. For 
HD 45088, OQ3 = H3 is dominated by OP3 = K 3 . In 
contrast, for £ TrA, OQ3 is dominated by P3Q3 = 63. 
This difference accounts for detection of K3 for HD 45088 
and its non-detection for ( TrA. 

Despite the enormous statistical significance of the e 
reported by Skuljan et al. (2004), it fails to pass the tests 
of reality developed in Sect. 3.2 and 3.3. Thus follow-up in- 
vestigations predicated on the reality of e carry substantial 
risk of misleading or erroneous conclusions. Nevertheless, 
just as p\ > 0.05 does not prove that an orbit is not eccen- 
tric but merely that e is not significant, so a finding that 
P2 < 0.05 and p 3 > 0.05 does not prove that e is not real 
but merely that this supposition does not find significant 
support in the data. 



8 



Lucy: Elliptical orbits 



6. Conclusion 



Appendix A: Fourier coefficients 



The aim of this paper has been to develop tests that allow 
SBs with elliptical orbits of high reliability to be sepa- 
rated from those with spurious e's. The example of HD 
45088 shows that the literature contains orbits that are 
not challenged by these tests. But £ TrA shows that some 
published e's are, as expected, probably due to perturba- 
tions. Evidently, extensive further recomputations would 
greatly clarify the reliability of the accumulated data on 
SBs. Indeed, categorizing SBs according to their location 
in the (pi,p2,P3)-cube offers an objective treatment of 
orbit quality that could replace or complement the sub- 
jective judgements (e.g., a, b, c, . . .) traditionally used by 
compilers of catalogues. Note that the recently published 
9th catalogue of SB orbits (Pourbaix et al. 2004) compiles 
individual velocities, thus greatly facilitating orbit recom- 
putations. 

As to measuring apsidal motion with high precision 
velocity curves (Sect. 1), observers should select SBs with 
Pi < 0.05 and P2 > 0.05, and then aim for a new orbit 
that satisfies the first and if possible also the second of 
the inequalities given in Eqs. (16) and (17). By doing so, 
tests for the detection of K3 and possibly also K4 can be 
carried out. Such detections would establish the reality of 
e and thus justify 2nd epoch observations. 

In contrast, re-observing systems with p\ > 0.05 in 
the hope of discovering hitherto undetected e's runs the 
following risks: that e remains undetectable; that e is sig- 
nificant but a perturbation is detected; that e is significant 
and no perturbation is detected but e is too small for de- 
tection of K 3 . Evidently, this (re-)observing strategy will 
have a much smaller discovery rate of orbits in the highest 
category of reliability (Sect. 3.4.1). 

A long standing problem for orbits of SBs is their non- 
uniform distribution of ui, sometimes referred to as the 
Barr effect (sec, e.g., Batten & Ovenden 1968). This was 
present in the orbits published by LS, both for those with 
and those without significant e's, and was attributed to 
perturbations. Accordingly, it will be of interest to re- 
examine the distribution of u for orbits with and without 
perturbations according to the p2-test. An obvious predic- 
tion is that orbits with P2 > 0.05 will have a negligible or 
much reduced Barr effect, and those with p 2 < 0.05 will 
have the effect enhanced. 

Given the typical precision of published velocity 
curves, a programme of recomputation would investigate 
the effects of perturbations with amplitudes <; 0.1km s _1 . 
But the work of Skuljan et al. (2004) demonstrates that 
perturbations are detectable to much lower amplitudes. 
In this regard, it would be interesting to observe eclips- 
ing binaries with light curves that imply ecosu; = to 
high precision. Comparison with the spectroscopic values 
of ecosu could then reveal, or demonstrate the absence of, 
perturbations down to ~ 1 m s _1 . 



The Keplcrian velocity of a binary component can be 
expanded in a Fourier series in the mean anomaly M 
(e.g., Plummer 1908). This is readily transformed into 
the Fourier series in the mean longitude L — M + lu 
given in Eq. (2). If we define a = K (1 — e 2 )cosiv and 
(3 = K\J\ — e 2 sinuj, the Fourier coefficients are 



and 



Cj = a Cj cos juj + (i Sj sinjui 



Sj = a Cj sin jui — (3 Sj cos jui 



(A.l) 



(A.2) 



Expressions for the coefficients Cj and Sj in terms of Bessel 
functions are given by Plummer (1908). Expansions of the 
latter as power series in e (e.g., Smart 1953, p. 378) yields 
the following formulae, accurate to 0(e 7 ) 



ci 



si 
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s 2 = e 
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where 



1 



Cj - 1)1 



(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.10) 

(ATI) 



Values of Cj and Sj obtained with these expansions 
have been checked against values obtained by numerical 
integration of the Fourier integrals. 

References 

Batten, A.H., & Ovenden, M.W. 1968, PASP, 80, 85 

Griffin, R.F., & Emerson, B. 1975, Observatory, 95, 23 

Lucy, L.B., 1989, Observatory, 109,100 

Lucy, L.B., & Sweeney, M.A., 1971, AJ, 76, 544 

Lucy, L.B. & Sweeney, M.A., 1973, Observatory, 93, 37 

Monet, D.G., 1979, ApJ, 234, 275 

Plummer, H.C., 1908, ApJ, 28, 212 

Pourbaix D., Tokovinin A.A., Batten A.H., Fekel F.C., 

Hartkopf W.I., Levato H., Morrell N.I., Torres G, Udry 

S., 2004, A&A, 424, 727 
Press W.H., Teukolsky S.A., Vetterling W.T., Flannery 

B.P. 1992, Numerical Recipes. Cambridge Univ. Press, 

Cambridge 



Lucy: Elliptical orbits 



9 



Russell, H.N. 1902, ApJ, 15, 252 

Skuljan, J., Ramm, D.J., & Hearnshaw, J.B. 2004, MNRAS, 
352, 975 

Smart, W.M., 1953, Celestial Mechanics (Longmans Green, 

London), p. 378 
Wilsing, J. 1893, A.N., 134, 90 



